Spin relaxation in a two-electron quantum dot 
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The electronic structure of artificial atoms in quantum dots (QD's) [f] can be directly 
studied in transport measurements. Such voltage-tunable experiments in a varying magnetic 
field (e.g., see Refs. [1-7] and the works cited therein) enable to change the number of 
electrons confined in a QD as well as the mutual arrangement of electronic energy levels. 
The arrangment of levels depends on the total spin of electronic configurations [3,4]. The 
application of the magnetic field induces transitions in the ground state between states 
characterised by different total spin. In this paper we focus on a role played by total spin 
in a simplest system, a two-electron droplet. At low magnetic fields the ground state is a 
singlet with total spin S = while at higher magnetic fields the ground state is a triplet 
with total spin S = 1, in analogy to the parahelium-orthohelium transition , already studied 
theoretically in Refs. [4,8] . Both singlet and triplet states and the singlet-triplet (ST) 
transition have been observed experimentally in both vertical [2,6] and lateral QD's [5,7]. 
However, while at least one of the triplet excited state is observed for magnetic fields below 
the ST transition, the singlet excited state is not seen past the ST transition, resulting in 
asymmetric (in B) high source and drain transport spectra. In the present work we present 
model calculations of total spin relaxation due to the mixture of spin-orbit (SO) and electron- 
phonon interaction which helps to explain the unusual behaviour of the levels associated with 
different total spin seen in transport experiments. 

We start with the SO interaction Hamiltonian for a two-dimensional (2D) electron in a 
quantum well written as [9-13]: 

H so = a(kx&) z -(3(k-&) , (1) 

where the layer plane is determined by the principal axes (x, y) of the crystal. This expression 
is a combination of the Rashba term [14] (with the coefficient a), and of a 2D analogue 
of the Dresselhaus term (with the coefficient (3) [10,11]. We use the following notations: 
k = — iSJ + eA/c and k = (k x ,—k y ) are 2D vectors; a x ' v are the Pauli matrices. The j3 
coefficient is determined by the formula [10,11]: 

(3 = %o^£== ■ Ry* • a 2 /d 2 = B-Rf- a 3 /d 2 , (2) 

where G is the band-gap width of the semiconductor, 7 so is the spin-orbit constant 
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[15], a = h 2 e/m*e 2 is the effective Bohr radius. (For GaAs: G = 1.52 eV, 7 so = 0.07, 
a = 9.95 nm Ry* w 5.74 meV, m* = 0.067m e , and B = 0.0043.) The parameter d is 
determined by averaging the square of the wave- vector ^-component of a 2D electron in the 
layer: d~ 2 = — (f \d 2 /dz 2 \ f) , where f(z) is the corresponding size-quantized function. (In 
the well we consider that / = ^/tt) 1 / 4 ^ 22 ^ 2 /y/d.) 

The Hamiltonian of the system of two interacting electrons in the harmonic potential 
^m*iu 2 (x 2 + y 2 ) may be written in Center-of-Mass (CM) and relative coordinates as H = 
H m (r) + CKm(R) +T~Ls- Here R = (ri + r 2 )/2 and r = ri — r 2 are Center-of-Mass and relative 
(Jacobi) coordinates of two particles. The first term in the Hamiltonian is 

(here cu c = eB/m*c; /i = m*/2 is the reduced mass, h = uj c /u is the dimensionless magnetic 
field). The expression for IKm(R) may be found from Eq. (3) for H m after the substitution: 
r — > R, \i — > M — 2m*, m — > M, and e — > oo. 

The term Tis is the spin-dependent part of the Hamiltonian, namely it is a combination 

For two 

electrons, the SO part + H s 2 q can also be written in CM and relative coordinates R 
and r. Denoting E = <xi + <x 2 , <r = <Xi — <r 2 . we obtain 

^5 = (ffflso + H rSO ) + ^ BJ BE 2 , (4) 



of the Zeeman and spin-orbit coupling terms: H,s = J2i=i,2 \9V»bB\oI + H$q 



where 



and 



Hrso = -D+ (at. + i/3E+) - D_ (at + - i/?E_) , 



H r so — —d+ (aa_ + i/3a + ) — <9_ (aa + — . (5) 

The new operators are 

[(X, Y) and (x, y) are the components of 2D vectors R and r] , b = m*u c /2h and 



The wavefunction of the two-electron system may be written in the form \l>Af(-R)e jM<I> • 
^ , m(^)e im<?f 'Xo-i,(T2 5 where a 1 and a 2 are the spin variables of the electrons. We expand the 
wavefunction in the basis set of the singlet and triplet states 

|s)=*oWo(r)|0,0) and |t, S g ) = * Wi(r)e^|l, S z ) (S Z = 0,±1). (6) 



Here fy is the ground state function (i.e. it obeys the equation IKo^o = h^oy 1 + h 2 /4& ), 
while the functions tpo,i have to be found from the equations 

H rn i) rn = E m i) m . (7) 

The analytical solution of Eqs. (7) could be found if I = ^Jh/m*uJo(h 2 + 4)~ 1//4 <C a 
or I ^> a . The / <C a case has been studied perturbatively in Ref. [8]. Here we con- 
sider the opposite limit I ^> a (this seems to be more relevant to a typical experimental 
situation). Then in the leading approximation the solutions of Eqs. (7) are the states of 



one-dimensional oscillator with mass ji and frequency ujq\Jz + 3/i 2 /4 localized in the vicin- 
ity of r = [2/g/a (l + /i 2 /4)] 1 ^ 3 3> Iq (here we have designated l = \Jh/mFuJo). In this 
approximation the energy measured from the ground state is 



h 2 m 2 h . . _ . Tiujq r 
5 m ,s z = E m -E = -g—; - - -u c [m + g S z ) = —D m ^ Sz (h) (8) 

TqTYI Z 4: 



(g* = gm*/m e ps 0.029). The equation determining the ST crossing takes thereby the form 
8i,Sz = a o/^o = 1/3; then for the ST crossing point we obtain h = 0.64 at S z = 0, 
whereas for the same ao/lo and S z the exact numerical calculation [8] gives h = 0.69). 

We now turn to the effect of SO interaction on the mixing of singlet and triplet states. 
Operators E± and S z commute with E 2 , therefore the first and third terms in Eq. (4) are 
not responsible for mixing of the singlet and triplet states. On the contrary, the second 
term in Eq. (4) results in this mixing. Indeed, let \S,S Z ) be the normalized spin states of 
two electrons. When operating on the state |0,0), the term H rS o yields the following non- 
zero matrix elements: (1, l\H rSO \0, 0) and (1, —l\H rS o\0, 0). (One can check that <r±|0, 0) = 
T2- 1 / 2 |l,±l).) 

Hence we see immediately the |0,0) singlet state is mixed with the |1,±1) triplet states 
but not with the |1,0) state. This state is therefore long-lived. 



For the states which are coupled, the expansion in terms of the small parameter oq/Iq 
leads to the following result for the mixing matrix element between \s) and \t, 1) states: 

M so = (s\H rSO \t, 1) « (I + 6r ) = ■ (9) 

If we take into account the Rashba term, we would find that another non-zero matrix element 
is (s\H r so\t, —1), however the state \t,0) is never mixed with the singlet. Further we neglect 
the Rashba coupling (usually there is a < (3 in GaAs hetero-structures, besides a vanishes 
in the case of an ideally symmetric quantum well). 

With help of Eq. (9) we find the hybridized states |S) = Cq\s) + Cf \t, 1), and |T) = 
Co\s) + where 



M. 



so 



2 M, 



sol 



1,^ 
+ SE 



1/2 



,± , ! M so 



Cf = ± 



2 M, 



sol 



SE 



1/2 



(10) 



(SE = \fsi~L + 4 |M 50 | 2 ^ • The corresponding energies are E T / S = E + (5i,i ± SE) /2. 

The next step in our study is the calculation of the |T) — > |S) (or |S) — > |T)) relaxation 
rate for the case, where S± t i > (or S^i < 0). Evidently the main relaxation channel is 
determined by emission of the acoustic phonon with energy hc s k = SE; where k = (q, k z ) is 
the phonon wave vector, and c s is the mean sound velocity (we use the so-called isotropic 
model, i.e., c s does not depend on the polarization and on the k direction; we consider that 
c s = 3.37 ■ 10 5 cm/s). The probability of this event is determined by 

■ = £ "^^ -6), (11) 

where M(q, k z ) is the appropriate matrix element 

M(q, k z ) = (T|C^ph|S> = C +*CV ((s|C^ph|s) - (l,t|f/ c -ph|t, 1)) • (12) 

Here the phonon field is averaged over the angle = rAq: 

V2 / fc \ 1/2 



(qr) 



£4- P h(R, r) = Us(d, fc,)e*» R (e** + e -*') = 2 (yj £ U s (q, fc,)e* R J 

(13) 

[s is the polarization, V is the sample volume, and £/ s (q, fc 2 ) is the renormalized vertex which 
includes the deformation and piezoelectric fields created by the phonon]. The integration 



with respect to z has already been performed, and reduces to the renormalization U s = 
U s {<\,k z )F{k z ), where the formfactor is F{k z ) = (f\e ihzZ \f). By using Eqs. (6), (12) and 
the expansion Jo(qr) m J(qr ) — q(r — r )Ji(qro), we obtain the matrix element (12); and 
after the substitution into Eq. (11) we find that the relaxation rate is proportional to the 
| Y, S U S \ 2 - The latter is represented as [12] 

nhc s k , 1 1 5pl 



Po t a 



where 



— + 



(14) 



r A (k) T D k & T P 

The summation involves averaging over the directions of the polarization unit vectors for 
both components of the electron-phonon interaction. The nominal times for the deformation 
and piezoelectric interactions in GaAs are r D ps 0.8 ps, and r P fa 35 ps [12,16]. The nominal 
momentum is p = 2.52 • 10 6 /cm [12,16]. We also refer to Refs. [12,16] for details concerning 
the meaning and the expressions of these values in terms of the GaAs material parameters. 

Finally, with the help of Eqs. (2),(9),(10),(12)-(14) and (11) we calculate the relaxation 
time. 

\/ T {h) = W-SC 2 - exp (-AS 2 ) ■ J . 

where 



7r(Ry)X°£ 2 



45 i 

£(h) = t— = ^Vlih) + 2 \BaMh)/d 2 } 



216(hc s ) 3 poTplod' 1 ' hu 
[see the definitions for C and 2\i in Eqs. (8-9)], 
'dalR Y \ 2 ^ /-i da 



A 



4llhc s 



' j m = L 

Jo 



5£-fe + S£ 2 (h) 



<V(h) 



n(O = a Ry/hc s (a /l ) 2/3 ^S(h)(4 + h 2 )- 1/3 , S = (t p /t d ) (a 2 Ry/2l 2 o nc sPo ) 2 , 



and 



V{h) 



al/llV^+h 2 - - (da /ll) 2 ] [S(h)a Ry} 2 /8(hc s ) 



As an illustration Fig. 1 shows the relaxation rate as a function of the magnetic field on 
the logarithmic scale (the main picture) and in the usual scale (the inset). The relaxation 
time is seen to have a sharp maximum in the vicinity of the ST crossing but constitutes 
a comparatively small value (of the order of 0.1 mcs) in the regions where the singlet and 
triplet lines are resolved. The non-monotonic behaviour of r on the right of the ST crossing 
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FIG. 1. The ST relaxation time r calculated for a /l = 1/3 and d = 5nm. The maximum 
corresponds to ST transition at u> c /ujq = 0.64. See the text for details. 

originates from the correlations between the wave-function characteristic distance r and the 
wavelength hc s /SE of the emitted phonon. 

We now turn to the discussion of the manner in which the ST relaxation could influence 
the transport spectroscopy through the QD states. By studying the kinetic processes of 
filling and emptying the dot in the presence of a large "source-drain" voltage, we estimate 
the effective electron life-time inside the dot at the "working level" i.e. at the level which 
participates in the transport through the dot. This effective life-time Tdot under experimental 
conditions of Refs. [6,7] we estimate to be of the order of 1 mcs, and this value should 
evidently be compared with the ST relaxation time calculated above. If the working level 
is exactly the upper level of the two-electron droplet, then the relaxation could influence 
the current. Namely, if r < r^t, then the working level could be emptied due to the ST 
relaxation occurring within the dot. In this case the current through the upper two-electron 
state becomes negligible. The relaxation process is asymmetric across the ST transition. 
Before the transition the "working level" involves \t, ±1) and \t,0) [see Eq. (6)] triplet 
states. (The Zeeman splitting is not resolved.) The \t, 0) state is long lived and hence 
observed in experiment while the \t, 1) state relaxes efficiently to the S = singlet state 
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\s). Past the ST transition, the \t, 1) state is the ground state but the excited state is the 
singlet. The singlet state relaxes efficiently to the \t, 1) ground state. Hence this asymmetry 
in the relaxation processes associated with the singlet-triplet transition could be responsible 
for the anomalies observed in transport experiments [6,7]. 

In closing it is worthy to mention other relaxation channels which are not taken into 
account in our calculation but which in the framework of the considered mechanism could 
additionally reduce the ST relaxation time. These are provided by special phonon modes 
(e.g., by surface and confinement phonons excited in the hetero-junction) and certainly by 
the SO Rashba coupling if the latter is significant. 
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